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Experimental evidence for superionic water ice 
using shock compression 
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Federica Coppari', Dayne E. Fratanduono', Raymond Jeanloz?, Damian C. Swift’ and Jon H. Eggert’ 


In stark contrast to common ice, l, water ice at planetary interior conditions has been predicted to become superionic with fast- 
diffusing (that is, liquid-like) hydrogen ions moving within a solid lattice of oxygen. Likely to constitute a large fraction of icy 
giant planets, this extraordinary phase has not been observed in the laboratory. Here, we report laser-driven shock-compres- 
sion experiments on water ice VII. Using time-resolved optical pyrometry and laser velocimetry measurements as well as sup- 
porting density functional theory-molecular dynamics (DFT-MD) simulations, we document the shock equation of state of H,O 
to unprecedented extreme conditions and unravel thermodynamic signatures showing that ice melts near 5,000 K at 190 GPa. 
Optical reflectivity and absorption measurements also demonstrate the low electronic conductivity of ice, which, combined 
with previous measurements of the total electrical conductivity under reverberating shock compression, provides experimental 


evidence for superionic conduction in water ice at planetary interior conditions, verifying a 30-year-old prediction. 


dinary polymorphism, with subtle rearrangements of the H,O 

molecules induced by pressure and/or temperature yielding at 
least 18 crystalline and several amorphous ice phases’*. Under high 
pressure — above 50-100 GPa — and near three-fold compression, 
water molecules lose their characteristic 100° H-O-H geometry 
and the cubic ice VII structure (or ice VIII at low temperature) 
transforms into the ionic ice X with linear O-H-O bonds™. Early 
classical molecular dynamics (MD) simulations on water ice under 
high pressure-temperature conditions? suggested the occurrence 
of superionic conduction enabled by liquid-like hydrogen ions 
rapidly diffusing within the solid lattice of oxygen. Several groups 
have since then reported Car—Parrinello and Kohn-Sham density- 
functional-theory-based MD simulations (DFT-MD) showing that 
the onset of superionic behaviour might be occurring at higher 
pressure®'°; these provided detailed predictions of its atomic struc- 
ture’ and transport properties'*'”"*. In particular, the high diffu- 
sivity and effective charge ~+e of superionic ice protons!’ enable 
high ionic conductivity o, exceeding 100S cm” whereas the elec- 
tronic transport remains that of a semiconductor, with electronic 
conductivity o, « o; Theory also predicts a rapid rise of the melting 
line associated with the onset of superionic conduction, reaching 
temperatures higher than previously thought”"'’, thereby suggest- 
ing that superionic ice could be stable under the expected extreme 
conditions inside Neptune and Uranus. 

In parallel with these computational efforts, static compression 
experiments at high pressures and temperatures” indicate a sud- 
den increase of the melting line above 20-50 GPa (refs **°’) with 
spectroscopic fingerprints” pointing to the transition to superionic 
ice. However, unambiguous thermodynamic and transport signa- 
tures of the superionic state remain elusive. Despite a claim of supe- 
rionic conduction”, electrical conductivity measurements up to 
900 K and 70 GPa also fall short by two orders of magnitude in o of 


| xtensive studies of the unique properties of water reveal extraor- 


identifying the expected characteristic high conductivity of rapidly 
diffusing protons. 

Dynamic compression offers an alternative approach for recreat- 
ing the extreme conditions of interest in the laboratory. However, 
single-shock compression of ambient-condition liquid H,O is 
not suitable to reach the predicted stability domain of super- 
ionic water because shock temperatures rise too quickly (Fig. la). 
Measurements of the optical properties reveal the transition, on 
increasing pressure, from transparent’? to opaque near 50GPa 
(refs °?) and finally to reflecting***° above ~100 GPa and 5,000 K. 
Given that ions do not contribute to the optical conductivity at such 
high photon frequencies ~10'*Hz, these optical measurements are 
probing the electronic contribution to the conductivity (o,) and the 
high shock-front reflectivity reveals the onset of significant elec- 
tronic conductivity in the dense fluid with o,>100S cm‘. In con- 
trast, electrical impedance measurements of the total conductivity o 
show the emergence of ionic conductivity’ in transparent, shock- 
compressed water, with a rapid increase from 0.1S cm"! at 10 GPa 
to 20S cm at 20 GPa before saturating near 30S cm"! above 30 GPa 
and 2,000K. This saturation has been interpreted as a signature of 
water being fully dissociated: increasing pressure and temperature 
do not create more charge carriers, in agreement with some” but 
not all” DFT-MD simulations. 

The electrical (total) conductivity o was also measured using a 
series of reverberating shocks between 70 and 180 GPa, in two inde- 
pendent studies“, In these experiments, dividing the dynamic 
compression into multiple stages reduces the production of 
entropy and keeps the temperature below ~5,000 K while exploring 
higher density states (Fig. 1a). These reverberating shock experi- 
ments reveal electrical conductivity o~150S cm”, much larger 
than o~30S cm" at similar temperature along the single-shock 
Hugoniot (locus of shock states). This indicates that higher pressure 
and density induce a significant increase of electrical conductivity, 
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Fig. 1 | Experimental approach and representative ultrafast line-imaging pyrometer (SOP) and interferometric Doppler velocimeter (VISAR) data. 
a, Simplified phase diagram for water shows the predicted stability domain of superionic ices and the pressure-temperature states accessible by single 
(black line) or reverberating*' (blue diamond) shock compression of liquid water with pọ ~1g cm? or single shocks on pre-compressed ice VII with 
Po~1.6g cm” (blue line; see also Supplementary Fig. 10). b, Sketch of the diamond anvil cell (DAC) to pre-compress water into ice VII at 2.5 GPa. 

c, Schematic detail of the pressure chamber. d, VISAR and SOP data for shot 70998 overlayed with the corresponding shock-front velocity history 
U.(t) (green) and averaged SOP counts lineout (red) showing: (1) shock arrival at the diamond/quartz interface; (2) propagation through the quartz; 
(3) arrival at the quartz/ice VII interface; (4) propagation through the ice; and (5) arrival at the ice/diamond interface. Event (3) is used as t=O. 


reflecting either a higher ionic conductivity or the onset of a signifi- 
cant electronic contribution. Such an increase of ionic conductivity 
seems unlikely in the fluid phase, considering that o saturates above 
30 GPa along the liquid Hugoniot. Instead, the efficient diffusion of 
H ions in superionic ice is expected to enable high ionic conductiv- 
ity o,>100S cm", thus it has been proposed that these reverberat- 
ing shock experiments validate the prediction for superionic ice'”"*. 
However, experimental data on the optical properties of H,O at 
these conditions are critical to understand whether the emergence 
of significant electronic conduction 6, > 100 S cm"! may offer a valid 
alternative interpretation and argue against the existence of super- 
ionic water ice. 

Here, we report time-resolved optical pyrometry and laser 
velocimetry measurements at the Omega Laser facility, as well 
as supporting DFT-MD simulations of H,O in and beyond the 
pressure-temperature (P-T) regime explored by the reverberat- 
ing shock experiments (Fig. 1). Using a combination of static pre- 
compression and laser-driven shock compression, we document 
the shock equation of state of H,O along the Hugoniot starting 
with ice VII and find pressure-density and temperature—pressure 
signatures for melting near 190GPa and 4,800K (Fig. 2a,b). We 
also document the optical properties through reflectivity measure- 
ments (Fig. 2c) at 532 nm (2.33 eV) and absorption measurements 
at 532nm and near 650 nm (1.9 eV) (Fig. 3) and find that the onset 
of significant electronic conductivity above 100S cm™ occurs in 
the fluid phase at higher P and T conditions than reached in the 
reverberating shock experiments (Fig. 4). Hereafter, we refer to the 
optical properties as being for an energy of 2eV when we are com- 
bining reflectivity and absorption results. Altogether, the new data 
provide experimental evidence for superionic conduction in water 


ice at planetary interior conditions (Fig. 5) in agreement with the 
DFT-MD simulations. 

Water is pre-compressed into the cubic ice VII phase in a dia- 
mond anvil cell at room temperature, then further compressed and 
heated by laser-generated shocks with pressures between 100 and 
300 GPa (Fig. 1a). In these experiments, the shock waves are strong 
enough to transform the quartz reference plates into an opti- 
cally reflecting state with a submicrometre optical depth so that 
the velocimetry interferometer system for any reflector (VISAR) 
tracks the shock front velocity (U;) and reflectivity (R) while the 
streaked optical pyrometer (SOP) measures the thermal emission 
of a thin layer of shock-compressed material just behind the shock 
front (Fig. 1 and Supplementary Fig. 2). However, we find that the 
optical depth of shocked H,O is larger, decreases with increasing 
shock pressure and only becomes submicrometre above ~200 GPa 
(Fig. 3). Nevertheless, even at low shock pressure, we can track 
the shock velocity by taking advantage of the 1-2% reflectivity 
caused by the sharp refractive-index jump between unshocked and 
shocked ice. 

Using the Rankine-Hugoniot conservation relations and imped- 
ance matching with the quartz as a reference, the velocity history 
extracted from the VISAR fringe phase analysis provides shock 
velocity-particle velocity, U;-u,, points, equivalent to a pressure- 
density shock equation of state (Fig. 2a and Supplementary Fig. 4). 
As expected, the shocks on ice VII explore higher density states 
than the liquid Hugoniot, and the compressibility is well captured 
by the DFT-MD simulations. In the absence of a phase transforma- 
tion, most materials exhibit near-linear, slowly varying U,—u, rela- 
tionships. Changes in volume and compressibility caused by phase 
transformations are often revealed by jumps and kinks in U,-u,. 
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Fig. 2 | Shock equation of state and reflectivity along the ice VII Hugoniot. 

a, Shock velocity versus particle velocity compressibility measurements 
obtained from the impedance-matching analysis, linear fit and piecewise- 
linear fit. Three different standard statistical analyses (Fisher's F-test, Akaike’s 
information criterion test and Bayesian information criterion test) comparing 
the results of these two fits indicate that the piecewise model is a better 
representation of the data, indicating a small discontinuity near U,~15-16 km s7 
in agreement with the DFT-MD model. b, A small discontinuity in the shock 
temperature as a function of shock velocity is also observed, providing a 
second thermodynamic signature for the melting transition near U; ~15km s7. 
c, The shock-front reflectivity stays below ~2% up to ~17 km s? (250 GPa), 
indicating that the refractive index n is similar to the index of ice at 

300 K (ref. °), and then increases quickly with the onset of significant 
electronic conductivity above ~17 km s7. Errors bars represent standard 
deviations at the 16 level and include random and systematic uncertainties. 


The new data along the shock Hugoniot for water ice VII (Fig. 
2a) are best fit by a piecewise model with a small discontinuity 
near U,~15—16km s~. This suggests the occurrence of a phase 
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Fig. 3 | Optical absorption of shocked water ice VII. a, Velocimetry 
(VISAR) and pyrometry (SOP) streaked images for shot 70997 (185 GPa) 
with corresponding shock-front velocity history U.(t) (green) and averaged 
SOP counts lineout (red), showing extra brightness after t=O due to the 
incomplete screening of the hot quartz thermal radiance by the shocked 
H,O. b, SOP signal as a function of the shocked water layer thickness 

Ax = fo |[Us(t’)-u, (t’)}dt. Once corrected by subtracting the shocked H,O 
self-emission contribution, the observed signal exhibits a characteristic 
Beer-Lambert exponential extinction. c, Absorption coefficient a and 
optical depth d=1/a of the shock-compressed H,O obtained from the 
pyrometry and a similar analysis of velocimetry data, prediction from 
DFT-MD and estimated values from the optical conductivity inferred 

from the reflectivity measurements. Horizontal dotted lines correspond to 
selected optical conductivity ranging from 1 to 1,000 S cm? (Supplementary 
Information). Absorbance increases with pressure, but the corresponding 
electronic optical conductivity remains well below 100 S cm” up to 

200 GPa. Errors bars represent standard deviations at the 1o level and 
include random and systematic uncertainties. 


transformation with a positive Clapeyron slope, consistent with 
melting and in agreement with the DFT-MD predictions. 

As the pre-compressed water ice VII is transparent, the pyrom- 
etry (SOP) data provide direct access to the shock temperature of an 
optically thick layer of shocked H,O, using the shock front in quartz 
as an in situ brightness reference and a grey-body assumption. 
In addition to measurements on arrival of the shock in the ice layer 
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Fig. 4 | Electrical conductivity of shock-compressed water along the liquid 
Hugoniot (black) and ice VII Hugoniot (blue). Experimental data obtained 
in this study are compared with previous electrical impedance*”“°, optical 
reflectivity***> and absorption*°*? measurements and simulations" along 
the liquid Hugoniot, and electrical impedance*°"' and simulations along 
the pre-compressed Hugoniot. Guide-for-the-eye lines indicate the the 
evolution of the total conductivity (solid) and electronic conductivity 
(dashed). The low electronic conduction observed along the ice VII 
Hugoniot shows that the difference in the total conductivity along the liquid 
Hugoniot and along the ice VII Hugoniot cannot be explained by the onset 
of significant electronic conductivity. Instead, superionic conduction is 
likely responsible for the high total conductivity along the ice VII Hugoniot 
below 5,000K. 


(chosen as t=0), we also collect continuous measurements by fol- 
lowing the slow decay of the shock waves during their travel through 
the H,O (Fig. 2b). Provided we include corrections to account for 
quantum nuclear vibrational effects (Supplementary Information), 
the simulated shock temperatures agree well with the experimen- 
tal data. In addition, the initial density increase gained by the pre- 
compression allows us to reach lower entropy states and to access 
the predicted pressure-temperature stability domain of superionic 
water (Fig. 5). 

In the absence of a phase transformation, materials respond to 
increasing shock velocity by a steady temperature rise to accommo- 
date the rise in internal energy (Fig. 2b). Latent heat and/or changes 
in the thermodynamic properties on a phase transformation are 
often revealed by discontinuities or plateau in the T—U, Hugoniot. 
Here, we identify the small discontinuity that we expect on the 
crossing of the melting line near U;~15km s~. Together with the 
U,-u, discontinuity, this indicates that water melts along the ice VII 
Hugoniot near U;=15km s~ and T= 4,800 K. Converting the shock 
velocity into pressure with the measured U,—u, relationship, we can 
plot our data in the T—P plane and find that the melting pressure 
is 190 GPa (Fig. 5). This also shows that the pre-compressed shock 
states match the pressure and temperature probed during the rever- 
berating-shock electrical conductivity experiments’! and that 
these probed the properties of ice (see Supplementary Information 
for details on our re-analysis of refs ®*“ P—T conditions using 
DFT-MD equation of state data). 

Since the VISAR directly tracks the shock, analysis of the fringe 
amplitude using the shock front in quartz as a reference provides a 
measurement of the shock-front reflectivity R (Fig. 2c). Knowing 
the index of the unshocked ice VII (n)=1.52), the behaviour of R 


documents the changes in the refractive index of shocked H,O 
induced by increasing pressure and temperature. Between 120 and 
240 GPa, the measured weak reflectivity indicates that the index 
is similar to the index of ice at 300K at the same density”. The 
rapid increase of the reflectivity (up to 11%) near 350 GPa suggests 
the onset of significant electronic conductivity o,>100S cm”, 
similar to the conductivity observed along the liquid Hugoniot. 
Assuming that the real part of the index continues to increase 
slowly with density we estimate o,(2eV) ~900Scm™ at 290 GPa. 
Post-processing the DFT—MD simulations to calculate the opti- 
cal properties of shocked H,O confirms this interpretation, and 
reproduces the observed reflectivity and the estimated optical 
conductivity: superionic ice is a wide-bandgap semiconductor 
while the fluid is electrically conducting (Supplementary Fig. 9). 
Both a Drude model and the DFT—MD simulations suggest that 
the d.c. conductivity is similar to the optical conductivity at 2eV 
(Supplementary Information). 

Further analysis of the pyrometry and velocimetry records also 
reveals that shocked H,O is partially transparent up to ~240 GPa. As 
shown in Fig. 3 and Supplementary Fig. 2, the drop of the pyrometry 
(SOP) signal at the shock breakout from the quartz into the ice VII 
exhibits a significant broadening due to the incomplete screening of 
the hot quartz thermal radiance by the partially transparent shocked 
H,O. A qualitative assessment of the SOP data (Supplementary Fig. 
3) indicates that the optical depth decreases with increasing pres- 
sure. Correcting the signal to account for the thermal emission of 
the shocked H,O, and using the velocimetry data to calculate the 
thickness of the shocked H,O layer at each time reveals a character- 
istic exponential decay behaviour (Fig. 3b), which can be fitted to 
extract the optical depth d or its inverse the absorption coefficient 
a (Supplementary Information). Including a mirror on the quartz/ 
ice interface for three shots allowed us to collect additional extinc- 
tion data from the velocimetry laser probe, yielding similar values 
of a in the same spectral range. This quantitative analysis confirms 
the qualitative trend of a gradual increase of the absorption 
coefficient from ~10°cm™ near 120GPa to ~2X10*cm™! near 
240 GPa (Fig. 3c). Here, the DFT—MD predictions underestimate 
the observed absorption coefficients, suggesting that defect-induced 
tails in the electronic density of states (which cannot be modelled 
accurately within the simulation framework used here) might be 
dominating the dielectric function near 2eV and masking the pre- 
dicted 4-6 eV bandgap of the superionic ice (Supplementary Fig. 9) 

Combining the results on the reflectivity (giving n(w)) and the 
absorption (giving a(@)) up to 240 GPa, yields model-independent 
optical conductivity at ~2eV as o(w) =€,ca(w)n(w) (where €, is the 
vacuum permittivity and c is the speed of light in vacuum). Figures 
3c and 4 show that o,(2 eV) « 100 S cm™ below 200 GPa and 5,000 K, 
and reaches 200 S cm” only above 240 GPa, well above the melting 
line. Comparing our data for o, together with the previous measure- 
ments shows that the difference in total conductivity at similar tem- 
peratures along the liquid Hugoniot (single shock’”-*’) and along the 
ice VII Hugoniot (very close to the reverberating shock states***') 
cannot be explained by the onset of significant electronic conduc- 
tivity. Instead, considering also the two thermodynamic signatures 
in U,—u, and in T—U,, indicating that water is likely solid at those 
conditions, we conclude that ice has a very high ionic conductivity 
between 2,000 and 5,000 K (100 and 190GPa) along the pre-com- 
pressed Hugoniot. Our study therefore provides experimental evi- 
dence for the existence of superionic water ice, almost three decades 
after the first report of fast protonic diffusion in pioneering molecu- 
lar dynamics simulations of water at extreme conditions’. 

Providing experimental evidence for superionic water ice, our 
study also confirms the rapid rise of the melting line above ~50 GPa 
predicted by the DFT—MD simulations and hinted at by several 
static compression studies**”’ (Supplementary Fig. 10). This rapid 
rise can be attributed to the enhanced entropy of the superionic 
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Fig. 5 | H,O phase diagram at planetary interior conditions.*° Measured temperature for shocked water ice VII and computed curve from DFT-MD 
simulations. Liquid water Hugoniot experiments”? and DFT-MD prediction. Proposed experimental (broad white line) and simulated melting 
lines’, solid phase boundaries and diamond anvil cell data (see also Supplementary Fig. 10). Calculated temperature for the electronic conductivity 
measurements using optical absorption and the total conductivity measurements using electrical impedance*°“! (Fig. 4). Three sketches illustrate the 
transitions from the insulating solid ice to the oxygen solid lattice and diffusing hydrogens of superionic ice, and then the conducting fluid. Errors bars 
represent standard deviations at the 1c level and include random and systematic uncertainties. 


ice. The signature of a phase transformation observed in ref. ” near 
1,400 K and 70 GPa is therefore likely due to the solid-to-superionic 
ice transition, not to melting. 

Considering the abundance of H and O in the Solar System, and 
the stability of H,O molecules, water is expected to be an impor- 
tant planetary building block. As previously suggested'*", but in 
contrast with the prevailing model of fully fluid interiors”, super- 
ionic ice could constitute a large fraction of the interior of Neptune 
and Uranus“. However, the presence of crystalline ice does not pre- 
clude large-scale convection within these planets, as we expect a low 
shear strength for superionic ice, based on simulations of the shear 
modulus of superionic ice’, or by analogy with the high plasticity 
of other crystalline superionic compounds”. Yet, to compensate for 
its relatively low conductivity, near 100 S cm”, vigorous convection 
with flow velocities similar to those stirring Earth’s liquid outer core 
(~10‘m yr") would be required to sustain a dynamo in a superionic 
ice mantle**"’. The internal stratification allowed by the high melt- 
ing temperature of the superionic ice might then favour a thin-shell 
dynamo in the ionic fluid upper mantle*’. This could help recon- 
cile the observed multipolar, non-axisymmetric magnetic field of 
these bodies with their expected composition and gravity field**°. 
In contrast, conducting fluid water could constitute a large fraction 
of sub-Neptune exoplanets** with surface temperatures exceed- 
ing 500-700K due to high stellar irradiance, such as GJ 1214b or 
Kepler 11c among many others, and contribute to generate much 
different magnetic field topologies typical of thick-shell dynamos* 
(Supplementary Fig. 19). 

Understanding whether water and water-rich C-H-O-N mix- 
tures are miscible or not and what are the properties of the stable 
phases at the relevant conditions appears as an obvious extension 
of this work. More generally, as astronomical observations and 
extreme-condition laboratory experiments improve in the near 
future, we are confident that tying planetary modelling*“*** to 
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the fundamental properties of planetary constituents should soon 
yield a better understanding of the diversity of icy planet formation, 
structure and evolution. 


Methods 

Methods, including statements of data availability and any asso- 
ciated accession codes and references, are available at https://doi. 
org/10.1038/s41567-017-0017-4. 
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Methods 

Pre-compressed samples. Double-distilled, deionized water was loaded in 
diamond anvil cells (DACs) with 200 x 200 um quartz plates either 15 or 30 pm 
thick and four ruby microspheres”, then compressed at room temperature to 

2.5 GPa to form ice VII. The liquid-ice VI and ice VI-ice VII phase transitions 
were monitored by visual inspections during compression. Pre-shot pressure was 
determined by ruby luminescence, allowing us to infer the initial density near 
1.6 g cm” and the refractive index near 1.52 based on previous studies 
White-light interferometry provided the pre-shot thicknesses of the quartz and ice 
layers. The thin diamond anvil was 350 um thick and had a 10-15 um polyimide 
plastic coating facing the laser drive on top of a 4 um Au layer, and a 100nm Ti 
layer on the sample interface. The 1,200-j»m-thick diamond anvil and the quartz 
plate had anti-reflection coatings. Some targets had an extra 100 nm Ti layer on 
half of the quartz/water interface. 


Laser-driven shock compression. The DAC target assemblies*“° containing the 
water ice were shock compressed by direct-drive laser ablation at the Omega Laser 
Facility, Laboratory for Laser Energetics (LLE) of the University of Rochester, NY 
(USA) using 700-2,000J of 351 nm UV laser in a 1 ns super-Gaussian flat-top 
temporally shaped pulse. Up to six laser beams incident at a 23° angle with the 
normal to the target were used (Supplementary Fig. 1). SG8-Flat distributed phase 
plates created a super-Gaussian spatial distribution matching the 900 m DAC 
opening“. The shock compression is destructive: each DAC is shot once. 


Ultrafast diagnostics. Our primary diagnostics were Omega’s velocity 
interferometer system for any reflector (VISAR) and streaked optical pyrometer 
(SOP). The VISAR offers a line-imaging, time-resolved record of the velocity of 
moving reflecting interfaces and optical properties (reflectivity, absorptivity) at 
the probe laser wavelength, 532 nm. The fringe phase shift encodes the velocity 
history while the fringe amplitude records the reflectivity of the tracked moving 
interface and the absorptivity of the different media that the probe laser is passing 
through®*. The SOP acquires line-imaging, time-resolved records of the spectral 
radiance between 590 nm and 850 nm (refs °""*'). Both diagnostics share the 
same optical telescope and image the same area of the target. VISAR sensitivities 
(vacuum velocity per fringe, vvpf) were 2.732 and 6.907 km s~. Two-gigahertz 
timing fiducials were recorded for each streaked image on all three cameras, 
allowing us to accurately correct the timing for sweep rate nonlinearities. The rise 
time of the SOP between 10% and 90% at a shock breakout event is ~130 ps. 


Equation of state simulations with DFT-MD. Using the Vienna Ab-initio 
Simulation Package (VASP), we performed finite-temperature DFT-MD 
simulations of water, extending previously reported work'®'*, to precisely map 
the pressure P and internal energy e as a function of density p and temperature 

T in the experimentally explored region. In addition to simulations over a coarse 
p-T grid approximately every 1,000 K (Supplementary Tables 5 and 6), we also 
performed a series of simulations every 100 K between 4,000 and 6,000 K at three 
densities close to the ice VII Hugoniot. To obtain meaningful results and avoid 
metastability issues, we analysed each MD snapshot and discriminated between 
the fluid and the superionic phases. We then fitted the temperature dependence of 
both internal energy and pressure along the three isochores for both the superionic 
and the fluid phases. We therefore obtained two grids of DFT equation of state 
(EOS) points for the fluid and superionic phases (Supplementary Tables 7 and 8). 
We applied a correction to the vibrational energies of the protons to account for 
their quantum nature®”’, and were able to solve the Rankine-Hugoniot equation 
e=e,+0.5(P+Py)(1/p)— 1/p) to obtain the predicted shock pressure, density 

and temperature starting with water ice VII. The calculation was performed at 

the Perdew, Burke and Ernzerhof (PBE) level of approximation to the exchange- 
correlation functional in the NVT ensemble (constant number, volume and 
temperature) with a plane-wave energy cutoff of 900 eV, a time step of 0.4 fs and 

a Nosé-Hoover thermostat using 162 atoms. Projector-augmented wave (PAW) 
pseudopotentials were used to account for the oxygen core electrons: we used 

the O and H in the VASP PBE PAW library with core radii of 1.52 and 1.10 

Bohr, respectively. The electronic density was constructed from single-particle 
wavefunctions by sampling at the T point of the Brillouin zone. The electron 
occupation numbers were taken to be a Fermi-Dirac distribution set at the average 
temperature of the ions. For the different densities and temperatures, we used a 
sufficient number of bands such that we had at least 40 bands with occupation 
numbers smaller than 10~°. The MD simulation was run for 20,000 steps with the 
last 16,000 steps used for averaging the internal energy and the pressure. 


Data analysis. Tracking reflecting shocks with VISAR to measure shock velocity 
and impedance match analysis to obtain pressure-density data. In the experiments 
reported here, the shock waves were strong enough to transform the quartz into 
a metallic-like, optically reflecting state so that VISAR tracked the shock-front 
velocity U; and reflectivity R while SOP measured the thermal emission of a thin 
submicrometre layer of shock-compressed and shock-heated material just behind 
the shock front". 

Although H,O was not always transformed into a metallic-like state with a 
submicrometre optical depth, the high dynamic range of reflectivity measurements 
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provided by the line-imaging VISAR at Omega allowed us to track the shock front 
as it propagated through the water ice VII layer (except for shot 71795). Even at 
low shock pressure, when the optical depth of shocked H,O was larger than several 
micrometres, the VISAR detected the 1-2% reflectivity due to the sharp refractive- 
index jump at the shock front between the unshocked and the shocked ice. 

Before extracting the phase map, some images were corrected by subtracting 
ghost fringes that are caused by spurious residual reflections at the ice/diamond 
or diamond/vacuum interface’. Having determined the time dependency of the 
shock velocity, U,(t) (Fig. 1 and Supplementary Fig. 2), we could extract Us quartz 
and Us waten the shock velocities right before and after the time (which we set as 
t=0) of the shock transit from the quartz into the water layer, using a linear fit 
over ~100-200 ps extrapolated to t=0. 

The Rankine-Hugoniot relations”, derived from conservation of mass, 
momentum and energy across a shock front, link the initial and final shock 
pressures, internal energies and densities with the shock and particle velocities U, 
and u,. For a given set of starting conditions Py—p,—é, all the shock end states lie 
on a curve in the pressure—density—energy space, the Hugoniot. 

Using the Rankine-Hugoniot relations, the measured change in shock velocity 
on the crossing of the interface between two materials — having different density 
and/or compressibility — can be used to infer the shock compressibility of the 
second material (sample) if the behaviour of the first material (reference) is known. 
Here we use quartz as a reference for the impedance-matching procedure™”. 

We explicitly include corrections to account for the pre-compression of the 
quartz reference”. 

For each shot, the impedance-matching procedure applied to the Us Quartz» Us water 
pair yields a U,-u, point, which is equivalent to a P-p datum, so that we obtained 
eleven U,-u, and P-p data points (Fig. 2a, and Supplementary Fig. 4 and 
Supplementary Table 1). We used a Monte Carlo routine to explicitly propagate 
random uncertainties associated with (i) the VISAR record, (ii) the data reduction 
to obtain the phase shift map and the initial pressure measurements, as well as (iii) 
systematic uncertainties associated with the VISAR instrument and (iv) the use of 
the quartz reference for the impedance-matching procedure. 

In the absence of a phase transformation, most materials exhibit a near-linear, 
slowly varying U,-u, relationship. Changes in volume and compressibility caused 
by phase transformations are often revealed by jumps and kinks in U,—u,. 

We show the new data along the shock Hugoniot for water ice VII 
in Supplementary Fig. 2a. We applied both a linear fit U;=a+ bu, 

(a=4.378 + 0.355km s~“, b= 1.362 + 0.0406); and a two-segment piecewise 

linear fit defined as U;=a, + b,u, for u, <c and U,=a,+b, u, for u,>c 

(a, =7.547 +1.16km s~”, b, =0.931+0.155, a,=5.762 + 0.443 km s”, b,=1.222++ 
0.046 with c=8km s`). Three different standard statistical analyses (Fisher's F-test, 
Akaike’s information criterion test and Bayesian information criterion test) 
comparing the results of these two fits indicate that the piecewise model is 

a better representation of the data. Analysis of the velocimetry data 

suggests the occurrence of a phase transformation with a positive 

Clapeyron slope. 

Note that the presence of the mirror coating on about half of the quartz/ice 
interface in shots 71000, 71001 and 71795 (Supplementary Fig. 1) caused artefacts 
in the VISAR data that prevented us from accurately measuring the shock velocity 
in the quartz. 


Reflectivity. While the VISAR fringe phase encodes the velocity information, the 
amplitude of the fringes contains information about the reflectivity of the moving 
reflector (here the shock front). We used the shock front during its travel across 

the quartz plate as an in situ reference based on previous absolute measurements at 
the Omega Laser*’”' but using corrections to account for the higher initial density 
of the precompressed quartz reference’. Doing so, we were able to determine the 
reflectivity time dependence R(t) and extract one R-P datum per shot for shocked 
H,O using an extrapolation of R(t) to t=0* and the pressure obtained from the 
impedance-matching analysis (Supplementary Fig. 2c). We explicitly corrected the 
fringe amplitude to compensate for small sweep-rate nonlinearities by normalizing 
it to a constant dwell time per pixel. We also performed a Hill-function fit to the 
shock-front reflectivity for water along the ice VII Hugoniot that we use to estimate 
the emissivity of shocked H,O that we need to infer shock temperatures from the 
SOP data: R= 1.35 +5.82U 8°" / (17.510?! + Us"? Another fit restricted to the 
data above U;=16.5km s™ was used for shot 86363: R = 69.51—10.40U, + 0.3808Us. 


Pyrometry. The spectral radiance space-time history recorded by the SOP allowed 
us to obtain shock temperatures provided one can attribute the observed radiance 
to a homogeneous part of the dynamically compressed target package. Note that, 
as for the reflectivity analysis, we explicitly corrected the SOP signal to compensate 
the small sweep-rate nonlinearities by normalizing to a constant dwell time 

per pixel. 

Quartz shock-compressed into the fluid phase has significant metallic-like 
electronic conductivity (optical depth is much smaller than 1 um) so that the SOP 
collects the spectral radiance of about a 1 optical-depth-thick layer of shocked SiO, 
at any given time during the shock transit in the quartz plate. We can then relate 
the SOP counts to the shock temperature, assuming that thermal emission is that 
of a grey body, with an emissivity e independent of wavelength, that is inferred 
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from the shock-front reflectivity: € = 1—R. We used corrections to account for the 
higher initial density of the pre-compressed quartz reference, compared to the 
original absolute measurements”. 

For Omega shot numbers 69347, 70421, 70424, 70998, 70999 and 86363, 
the SOP images show a sharp (limited by our temporal resolution) drop in 
radiance on shock breakout of the quartz into the water ice layer, followed 
by an approximatively linear decay up to the arrival of the shock at the ice/ 
diamond interface (Fig. 1 and Supplementary Fig. 2a). Here, the interpretation is 
straightforward and we can attribute the radiance between these two events to the 
thermal self-emission of the shock-compressed H,O. Using the fit of the measured 
reflectivity to calculate the emissivity of shocked H,O, we obtain the temperature 
time dependence T(t) for each shot from which we can extract a T-P datum using 
a linear extrapolation of T(t) to just after t=0 and the pressure obtained from the 
impedance-matching analysis (shown as blue circles in Figs. 2 and 5). In addition, 
we can also combine U,(t) and T(t) during the transit of the shock front through 
the water layer to obtain a continuous T(U,) record. As the shock decays slowly 
during its transit, we collect a series of shock temperature data along the ice VII 
Hugoniot in each shot (shown as shaded blue areas in Figs. 2 and 5). We emphasize 
that these decaying-shock data differ from shock-and-release measurements 
that would be obtained if one were to look at an interface between a hot, opaque 
shocked material and a less compressible transparent ‘window’ or the vacuum 
(see also ref. °). 

The SOP records for the shots at lower pressure (Fig. 3 and Supplementary 
Figs 2 and 3) are slightly more complex, as the SOP did not record only the 
thermal emission of the shocked H,O during the shock transit in the water ice 
layer. Instead, the drop of the pyrometry (SOP) signal at the shock breakout from 
the quartz into the ice VII exhibits a significant broadening due to the incomplete 
screening of the hot quartz thermal radiance by the partially transparent shocked 
H,O. A detailed analysis of these records provides a way to accurately determine 
the optical depth of the shock-compressed H,O as well as the shock temperature 
for a few shots in which we could create an optically thick shocked H,O layer. 

Using the velocimetry records we obtained U,(t) from which we can determine 
the shocked water layer thickness at each time Ax(t) = S [U,(t')—u,(t')] dt" 

(the shock front moves at U; but the quartz/ice interface is also put into motion, 
to velocity u,, by the shock). We could then transform the measured SOP time 
history, counts(f), into a record as a function of optical transit depth, counts(Ax). 
As shown in Fig. 3, such traces exhibit an exponential-like decay, but to obtain a 
more accurate determination of the optical opacity, we corrected these records 
by subtracting the estimated contribution of the shocked H,O self-emission. To 
do this, we first calculated the time dependence of the H,O shock temperature 
using the measured U,(t) together with the T(Us) relationship predicted by the 
DFT-MD, then used the quartz calibration for shock temperature to transform 
these estimated temperatures into SOP counts. The good agreement between the 
measured and predicted shock temperatures provides a good consistency check. 
The corrected SOP traces exhibit a clear Beer-Lambert behaviour and exponential 
fits from t=0 yielded the optical depth d=1/a, where a is the absorption 
coefficient. The corresponding pressure was obtained by interrogating the U,(t) 
record at t=0* and at the instant when the corrected SOP signal reached a steady 
state. The corresponding temperature was interpolated from the predicted T(Us) 
below 200 GPa and the averaged experimental data above 200 GPa. 

This analysis relies implicitly on a steady thermal emission from the quartz that 
has isentropically released from its initial shock state due to the shock impedance 
mismatch with the ice. One-dimensional Lagrangian hydrodynamic simulations of 
the experiments using the commercial code Hyades” indicate that the temperature 
of the hot quartz at the quartz/ice interface stays constant during the few 
nanosecond duration of the optical absorption measurements. We also made sure 
to avoid attributing the drop in quartz temperature to absorption by the water layer 
by calculating the SOP counts corresponding to the release temperature, estimated 
from the impedance-match analysis at t=0. 

For the few shots in which we could create an optically thick shocked H,O layer 
with Ax > 7—10d (70997, 70423 and 71000), we also extracted the portion of the 
record after the instant when the corrected SOP signal reached a steady state to 
obtain continuous T(U,) data. Here, we made sure to account for the smearing of 
the decaying-shock data by the approximately micrometre measured optical depth, 
but found that this effect was negligible due to the slow velocity decay. 

The targets for shots 71000, 71001 and 71795 featured an additional titanium 
mirror coating on approximately half of the quartz/water interface. Orienting 
appropriately the slit of the streak camera that represents the imaging direction 
y, we recorded VISAR and SOP for both the coated and uncoated region for each 
shot (Supplementary Fig. 1). Supplementary Figure 5 describes the experimental 
configuration and data interpretation for shots 71000 and 71001 near 180 GPa for 
which the shocked H,O layer optical depth was about 2 um. In the portion of the 
VISAR image corresponding to the portion of the interface with a mirror, two 
reflectors are tracked simultaneously: (i) a weak ~1-2% reflection at the shock 
front due to the mismatch of the real part of the refractive index between the 
unshocked ice VII and the shocked H,O; and (ii) a strong 37% reflection at the 
titanium-coated quartz/ice interface that is being attenuated by the two passes 
through the partially absorbing layer of shocked H,O. These two contributions 
are visible in the portion of the VISAR streaked image shown in Supplementary 


Fig. 5. During the shock transit through the diamond and the quartz plate, the 
VISAR signal reflects off the Ti mirror coating and the fringes are stationary 
(no motion). On event (3), the arrival of the shock at the quartz/ice interface, the 
fringe intensity drops rapidly while fringes associated with the weakly reflecting 
shock front become visible and appear slightly slanted, indicating a slightly 
decaying shock. 

Analysis of the temporal history of the fringe amplitude for the two VISAR 
channels indicates that this decay occurs on a longer timescale than the etalon 
delay for each VISAR. When plotted against the thickness of the shocked H,O 
layer, the rapid decay is found to be in good agreement with the expected Beer- 
Lambert exponential behaviour. The optical depth values are very close to 1.5 pm 
(Fig. 3). Note that the extinction appears faster on the VISAR images because the 
attenuation occurs twice (before and after reflection from the Ti mirror coating). 
Here, we use the nominal 37% reflectivity of the titanium layer at the quartz/ice 
VII interface as a reference to measure the reflectivity of the shock front. 

Supplementary Figure 6 describes the corresponding pyrometry data. Taking 
advantage of the line-imaging capability, we recorded SOP data for both the 
coated and uncoated portions of the targets. In both regions, two contributions 
are collected simultaneously by the SOP as the shocked H,O is not optically 
thick: (i) a weak self-emission from the shocked H,O; and (ii) a stronger self- 
emission by the shocked-and-released quartz that is being attenuated by the 
one pass through the partially absorbing layer of shocked H,O. Looking at the 
region with the mirror, we see no self-emission until event (3), when the thin 
Ti layer is quickly heated by the arrival of the shock wave in the quartz plate 
underneath, and the pyrometry signal in this portion of the record quickly 
matches the one in the uncoated quartz portion. Here, we could not accurately 
measure the shock velocity in the quartz because the fringe quality was not 
sufficient. Instead, we calculated it from extrapolating the measured shock 
velocity in shocked H,O to t=0 as well as a linear fit of the Us guartz—Us water data 
obtained from the targets without a coating (using only the four points below 
the melting line). Having determined Us quar» We used the SOP during event (2) 
to link the SOP signal to the grey-body temperature. A portion of the record 
corresponding to the late time of shot 71000 is reported in Figs. 2b and 5. We 
also analysed the SOP decay to obtain the absorption coefficient as describe 
above (Fig. 3). 

In this shot, the shocked H,O remains optically thin: the optical depth is 
comparable to the 34 um thickness of the sample layer. The VISAR record is 
dominated by the strong reflection at the Ti-coated quartz/ice interface, only 
slightly attenuated by the shocked H,O. The weak signal associated with the ~1% 
shock-front reflectivity is not distinguishable. Identifying events (3) and (5) yields 
the shock transit time through the ice layer, which allows us to obtain an average 
shock velocity U;=12.4+0.8km s~. Assuming the compressibility (we used 
the DFT-MD model from French'’), we can calculate the true particle velocity 
u,=5.89 + 0.16km s` that corresponds to a pressure at event (3) near 120 GPa. 
Analysing the fringe shift between events (3) and (5) yields the apparent interface 
velocity u,,,=6.56 + 0.16 km s~. Using also the unshocked refractive index 
n= 1.52, we can then calculate the real part of the refractive index of the shocked 
H,O at these conditions using” n=(u,,,—yUs)/(u, — Us) = 1.91 + 0.19. 

This value is in excellent agreement with the measured values at the same 
density at room temperature”. We also analysed the VISAR amplitude to obtain 
the absorption coefficient (Fig. 3). For this shot, the SOP record (Supplementary 
Fig. 3) shows a clear drop at event (3) consistent with the expected drop from the 
incident shock state temperature to the released state temperature. Analysing the 
(very small) decay of the SOP signal due to the absorption by the shocked H,O 
layer, we can also obtain a similar absorption coefficient. 


Data availability. The data that support the plots within this paper and other 
findings of this study are available from the corresponding author upon 
reasonable request. 
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